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I. INTB0EUC1I0H 


Ice is an insidious enemy. It attacks cn two flanks, 
adding weight to the airplane and at the same time 
ruining its aerodynamic shape, as ice accumulates, 
more and more power is necessary to maintain speed and 
altitude, and the pilot gradually finds himself forced 
to sacrifice first his speed and then tit by bat his 
altitude in a desperate struggle to stay airborne [ 1 ]. 


ibe effects cf ice growth cn the performance of an 
aircraft are felt primarily through aerodynamic penalties; 
a drastic reduction in and a qTATT and an increase m 

drag. The two classes of ice accretions, known as rime and 
glaze, are formed under different flight conditions. Biae 
ice is formed at low air temperatures and low velocities. 
The droplets freeze on impact and usually are found in 
flight through clouds with low liguid water content. 

Figure 1 shows an example of a rime ice accretion with its 
characteristic streamlined leading edge. 

Glaze ice on the other hand is formed at temperatures 
near the freezing mark and higher velocities, kith this 
type of ice growth, a phenomenon known as runback occurs. 
Rather than freezing on impact, the water droplets travel 
a short distance before freezing. The resulting shapes are 
of the type shown in Figure 1, with the characteristic 
horns. It is with glaze ice accretions that the greatest 


l 



aerodynamic losses are found, and it is this type of ice 
that is the subject of this paper. 

Classically, most of the work done in the study of 
icing, particularly that dene fcy the HftCA in the 1950's, 
was concerned with mechanical means of preventing or 
removing the ice, known as anti- or de-icing. Bcwever, 
with the increase in general aviation aircraft with 
smaller powerplants and lighter weight, a mechanical 
method of solving the ice problem is no longer acceptable. 
Retro-fitting the aircraft components with pneumatic toots 
or heating elements tend to increase the aircraft's 
weight, cost, and complexity. 

A letter approach would be to design the component 
itself with characteristics that would reduce the chances 
cf ice growth and the detrimental effects if growth does 
occur. This process has been investigated by Bragg [2,3,4] 
for rime accretions but no attempt has been made fer glaze 
ice conditions. Gray [5,6] derived an empirical formula 
for predicting iced airfoil performance degradation tut 
the correlation has not teen found to fit recent 
experimental data very well. 

Rhen trying to develop a method for evaluating the 
glaze ice problem, two phases must he examined. The first, 
a thermodynamic problem, deals with the prediction of the 
actual geometry cf the ice shape. The second is to 
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determine a scheme for analyzing the performance losses 
incurred once the geometry cf the ice has teen determined. 
The study described in this paper applies current 
potential flew methods to this problem. The approach 
discussed is not a final solution to the problem. Bather, 
it is intended as a first step in developing a glaze ice 
analysis method. Further investigation intc the properties 
of the flowfield in the region of the ice accretion is 
required before a complete scheme can be formulated. 
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II. SOHVEI CF IITEBA1UEE 


Host of the early investigations into tbe icing 
phencnenon were concerned with de-icing. The first cf 
these efforts was the development of inflatable de-icing 
boots by the B.F. Goodrich Coipany in the 1930*s. This 
concept is still in wide nse today. Befinecents hare 
reduced the boot in its deflated fora to the point that 
its presence barely affects the geoaetry of the wing. 

The first aajoi investigation into the icing 
characteristics cf various airfoils and the resultant 
aerodynaoic penalties was perforaed by the NACA in the 
1950*s [5]. Information was gathered on the 65A004, 

63A009, 0011, 65-212, and 63-015 airfoils. However, few 
correlations were drawn between the aerodynaaic penalties 
incurred and tbe shape and location of the ice accretion. 

The first aajor effort to draw these correlations was 
by Vernon Gray £5,6] in the 1960*s at the lewis Research 
Center. Gray developed an empirical eguation which relates 
known icing conditions with change in drag coefficient. 

The major testing was performed in the BASA lewis 6 * x 9* 
Icing Eesearch Tunnel on the HACA 65A004 airfoil. A wide 
range cf parameters were examined, including icing time. 


4 



airspeed, freestream temperature, liquid water content, 
cloud droplet impingement efficiency, angle of attack and 
leading edge radius of curvature. However, the correlation 
he developed from this study does not readily predict 
changes in lift coefficient and moment coefficient. An 
interesting facet cf Gray's correlation is the ability to 
mathematically grew the ice at a given angle of attack and 
then study the performance changes at another angle. 

Recent data however has shown that even though Gray's 
correlation reasonably predicts ACd at the angle the ice 
is grown, its accuracy drops significantly when the 
calculation is performed at another angle cf attack. 

Some interesting observations by laschka and Jesse [7] 
came from ether investigations in the Lewis Icing funnel. 
They observed that as the angle of attack is varied, many 
different ice shapes will be obtained. Also they noted 
that when the time of the icing encounter, t, is varied, 
the ice height will be approximately proportional to the 
value cf t, while the impingement limits are time 
independent. 

In order to begin quantifying the performance 
degradation due to ice, a scheme had to be developed which 
could predict the flowfield about the irregular ice shape. 
In 1568, Dvorak published a method to predict the 
development of turbulent boundary layers over rough 


5 



surfaces [8]. This approach is incorporated in his 
program, which was investigated in this paper. 

In addition to the roughness effects associated with 
icing, is the existence of a large separation tuhtle in 
the area of the ice shape. little research has teen done 
studying these laminar separation babbles. Host computer 
programs, such as the Eppler code [9], when they predict 
laminar separation, consider this simply a transition 
point between laminar and turbulent boundary layers. 
However Venkateswarl and flarsden [10] investigated laminar 
separation bubbles that occur at 60-70* chord. They 
developed a correlation to predict the size and shape of 
the laminar bubble. Also in 1976, Crimi and Beeves [11] 
studied leading edge laminar separation bubbles and 
developed a scheme to predict the onset of transition in 
the shear layer. 

In the late 70's and the present, icing research has 
increased with the work of Ingelman-Sundberg, Shaw, firagg, 
Gregcrek and others. Ingelman-Sundberg and Trunov [12] 
published a joint report from the Swedish-Soviet Rocking 
Group cn Scientific-Technical cooperation in the Field of 
Flight Safety. Flight test and icing wind tunnel studies 
were performed and the concept of simulated ice was 
developed as a means of investigating the aerodynamic 
effects of ice growths. 
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Shan [13], Eragg and Gregcrek [2,3,4] continued 
investigation in tie lewis Icing Eesearch Tunnel in the 
1980’s. Extensive data were collected on the lift and drag 
penalties of ice growths. Bine and glaze ice accretions 
were ncdelled using mahogany and pressure tapped so 
detailed aerodynamic data could he collected. Ibis work 
serves as the primary database for the analytical effort 
to be presented in this paper. 

Of particular importance to the glaze ice analysis was 
the work of Pfeiffer and Zuowalt [14] and Hclachlan and 
Karancheti [15] who investigated the flowfield around 
airfoils with highly deflected spoilers. Pfeiffer and 
Zuowalt utilized a splitter plate arrangement to visualize 
the separated zones created by the spoiler. 

lastly, Bristow [16] has developed an inviscid 
computer program which allows for input of mixed 
analysis/design boundary conditions. Eor example, the 
input to the pregram can consist of an airfoil with its 
geometry partially defined and a desired pressure 
distribution in the undefined region. The program will 
then held the input geometry fixed and design the 
remaining portion of the airfoil based on the input 
pressures. This program was particularly useful in the 
author’s investigation of the separation zone associated 
with glaze ice. 
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Ibis review of literature should give the reader a 
clear picture of the deficiency of direct investigations 
into the glaze ice problea. It is hoped that the study 
reported here will spawn continued efforts in this area. 


8 



III. POTENTIAI FLOS THEOBY AND PAHELLIBG HETHCDS 


la order to analyze the performance degradation that 
occurs due tc glaze ice accretions, a method for 
predicting the flcwfield aad therefore the pressure 
distribution of the iced airfoil must te developed. As a 
first step in accomplishing this task, current potential 
flow computer pregrams were investigated. These potential 
flow solutions are based on an incompressible, inviscid, 
and irrotaticnal fluid, for which the classical Bavier- 
Stokes Equation can be reduced to the Laplace Equation, 

V 2 cJ> = 0 (1) 

One scheme presently in use to solve this equation 
involves the distribution of surface singularities on a 
closed polygon which approximates the airfoil contour. 

This method is known as panelling. Examples of computer 
programs using this technique are; 1) Smetana, by F. 
Smetana, D. Sumney, N. Smith, and fi. Carder [17]; 2) 
Eppler, by B. Eppler and D. Somers [9]; 3) Dvorak, by F. 

A. Dvorak and F. A. Woodward [€]; and 4) Bristow, by D. B. 
Bristow £16]. The potential flew method of each of these 
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progtais will be discussed in this chapter 


Smetana 

The Smetana program approximates the airfoil geometry 
by a closed polygon* Vertices are placed cn the perimeter 
of the polygon (Figure 2). The velocity potential for each 
of these vortices can be expressed by: 


<f> = T tan" 1 Y-Yo (2) 

27 X-Xo 


where r is the vortex strength and (Xo,Yo) is the location 
of the center of the vertex. Ihis potential satisfies 
Laplace’s Eguation, which is linear and therefore the sum 
of any number of these potentials also will be a solution. 
The corresponding velocity expressions can be obtained by 
differentiation of the potential: 

u = a* = -r Y-Yo 

3X 27 (X-Xo) ^+(Y-Yo) 2 

(3) 

v = 3<t> = r X-Xo 

3Y 2tt (X-Xo) *+ (Y-Yo) 2 

The contributions of each vortex to the net velocity at a 
point (X, Y) can then be treated separately and summed. 
Therefore, the net velocity components, 7 and v, are: 
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V 


K 

= 1 \ (X-Xon) 

Ztt L (X-Xo{Jr+(Y-Yo N ) 2 


where K is the total nuaber of vortices and (Xo,Yo) is the 
location of the center of the Bth vortex. 

The boundary condition that aust be satisfied is that 
the flew aust be parallel to the airfoil surface, adding 
the contribution of the freestreaa velocity to the 
velocity components induced by the vortices, this 
condition can be written: 


v = ( d£\ 
U ro +u V^dX/wing 


tan a 


(5) 


If we denote the right side of this equation as and 
define 




(Ym-Yom) 

(X M -Xo N )4(Y I , r Yo N ) 2 


and 


( 6 ) 



_ (Xm-Xon) 

(X m -Xo n ) v +(Ym-Yo n ) 2 


we can write Eguation (5) as. 



B,= bliri+bl2r2+ ... +bDfrie 

ZTrU oo +anr 1 +ai2r2+ • • • +aiK r K 

i 

B T r= bKiri+bR2r2+ ••• +bRKrK 

27rU oo +a Kl r l +a K2 r 2 + • • • +a KK r K 

This set of eguaticns is then solved for the needed values 
of the vortex strength, r • 

The influence coefficients a MN and b^ N are solved for 
convenience at the midpoints of each panel. However, from 
the geometry, only K-1 values of the coefficients can be 
calculated unless the polygon is closed. The trailing edge 
point is then given two indices, 8=1 and 8=K. Then the 
system is determinant and can be easily solved. 

To satisfy the Kutta condition at the trailing edge, 
Smetana chose 


r i " - r K 

which still satisfies the reguirement that the circulation 
at the trailing edge is zero, since the trailing edge was 
denoted by the indices 8=1 and 8=1, the net vortex 
strength at the trailing edge is r.,+r K =0- Thus Eguation 
(7) contains K-1 distinct values of and K-1 values of B^ 
and is therefore solvable. 

Lastly, in order to obtain surface pressures, Smetana 
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uses the equation: 


P = P - lp(u 2 4-v 2 +2u.U on ) 

oo ' oo ' 

which is derived from the Bernoulli equation: 

P T = Pg + 1 pu 2 


(9) 


( 10 ) 


where u is the total fluid velocity and is calculated 
usinq the vector magnitude foraula: 


u = /(U ro +u) 2 +v 2 


(ID 


Eppler 

The Eppler program is very similar in construction to 
the Smetana code in that both utilize vertices to provide 
circulation and both satisfy the saoe flow tangency 
boundary condition. However, the Eppler code satisfies it 
on the actual input geometry points. Also, rather than 
applying a point vertex, Eppler distributes the vortices 
parabclically along each airfoil panel. The geometry of 
the panels is determined by a cubic spline fit of the 
input coordinates. The vortex strengths at the endpoints 
of each panel are solved for in the same manner as 
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Smetana 


The vorticity distribution between the panel endpoints 
is obtained from the equation: 

r(e) = e (l - (12) 

a \ ij v 

where % is the length of the panel, e is the local panel 
abscissa, and Yp is a parabolic vorticity factor. Ibis 
factor is calculated using the vortex strengths at the 
endpoints of the two surrounding panels, integration of 
the vortex distribution is then reguired tc evaluate the 
velocity contributions of each panel. 

Ihe Kutta condition is satisfied as in the Smetana 
program. Ihe requirement again is egual velocities on both 
sides cf the trailing edge and zero normal velocity with 
respect to the trailing edge bisector angle. Ihus, enough 
circulation is generated that the trailing edge becomes 
the rear stagnation point. 

Dvorak 


The airfoil contour is again represented by an 
inscribed polygon. However each pair of adjacent panels 
has a triangular distribution cf vorticity across it. The 
airfoil is thus modelled by a series of overlapping 
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triangular vertex distributions. It the leading edge of 
the airfoil, the strength of the upper and loner surface 
vortices are set egual to insure snooth flow. 

Ihe Kutta condition is satisfied by setting the 
strengths of the vertices on the trailing edge panels 
egual to zero. However, doing this reduces the system cf 
equations to be solved to H equations with N-1 unknowns. 

An additional unknown is added by applying a constant 
source distribution on the inside cf the airfoil surface. 
It should be pointed out that like the vortex strength of 
the trailing edge point used by Smetana and Eppler, this 
unknown source strength is always very nearly zero for 
airfoils with closed trailing edges. 

Bristow 

Ihe Bristow code is similar in design to Smetana and 
Eppler, however the singularities used on each panel are 
linear source and vortex distributions associated with the 
classical third identity of Green. One of the particular 
advantages of this method is believed to be its high 
numerical stability when used in the design mode of 
operation. 

Ihe vortex distribution generated is linear on each 
panel and the source distribution can be either piecewise 



constant or linear. This choice however shews little 
effect on the resalts obtained. The source strength at 
panel midpoints,^ , is found sinply from the following 
eguation: 


a. = V N-i + U co sin ( 0 i- a > (13) 

X V 


where Vm. is the prescribed normal velocity at the panel 

iN ip 

midpoints and 0^ is the local slope. With the prescribed 
source strength evaluated, it is left only to determine 
the total potential at a panel midpoint induced by the 
simultaneous action of the vortex and source 
distributions. 

The Bristow code has a unigue feature. It can perform 
mixed analysis-design problems. The user inputs fixed 
geometry regions and the desired surface velocities in the 
design region. The prograw iawediately satisfies surface 
continuity by stretching the input starting geowetry in 
the design region. Then an analysis only solution is 
obtained from the combined source-vortex singularity 
scheme. The geometry of the design region is then modified 
using a first order inverse method to minimize the 
difference between calculated and input values of 
tangential and normal velocities. This process is repeated 
until the convergence criterion is met. 
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Laminar Separation Bubbles 


One cannot deal with glaze ice accretions solely using 
potential flow methods. This is due to the presence of a 
laminar separation bubble which forms behind the glaze ice 
horn. Short laminar separation bubbles have very little 
effect on the integrated aerodynamic loads and most 
computer analysis programs assume that the bubble simply 
represents a transition point from laminar to turbulent 
boundary layers, however, laminar babbles of the type seen 
with ice accretions are sufficiently large that their 
effect cannot be neglected. 

Host of the work done on laminar separation babbles 
has been of an experimental nature. This is due to the 
difficulty in analyzing the interaction between the viscid 
and inviscid flow in the reverse flow region inside the 
bubble. In addition, evaluation of the transition point 
from laminar to turbulent flow in the free shear layer 
becomes more complicated. 

A diagram of a typical flow pattern observed with a 
separation bubble is shown in figure 3. The laminar 
boundary layer first separates from the surface yielding 
the region of reverse flow. Transition to turbulent flow 
occurs in the separated shear layer shortly before 
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reattachment. The region is divided by the streamline 
which separates from the surface and reattaches 
downstream, The area below the separation streamline is 
known as the recirculation region or separation bobble. 

Host of the examinations into the problem of laminar 
separation bubbles have used the classical boundary layer 
assumption that 

3P = 0 (14) 

3Y 

From Schlichting [18] however, it is noted that this term 
is of the order cf the boundary layer thickness. For most 
cases, this would be a valid assumption. However, the 
bubble behind a glaze ice accretion is much thicker than 
normal boundary layers and therefore the assumption that 
this term can be neglected may not be valid. In Chapter 5, 
the order of magnitude of this term is investigated. 

In an effort to analyze these separation bubbles using 
potential flow schemes, the assumption that the pressure 
gradient term is negligible will be considered valid. 
Pressures measured experimentally at the airfoil surface 
will be input to the Bristow code in the design mode. The 
corresponding calculated bubble shape will then be 
compared with the flow visualization results. 
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IV. EXPEBIHENTAL PBOGBAH 


Very little experimental data has been available on 
the performance degradation of airfoil sections resulting 
from ice accretions. To help alleviate this, a two-year 
test program was conducted in the HASA Levis 6* x 9* Icing 
Research Tunnel (IBT-Figure 4) . Its primary objectives 
were: 


1) . To examine a method of simulating ice accretions 
with wood shapes which were instrumented with surface 
pressure taps to obtain aerodynamic data. 

2) . To study and document the complex f lowfield in the 
region cf the ice shape through pressure 
distributions and flow visualization techniques. 

3) . To expand the current database of performance data 
on airfoils under icing conditions [4,19]. 


The first tunnel entry in 1981 was an actual ice 
accretion study. Glaze and Bime ice shapes were grown on a 
1.36 m chord NACA 63A415 model. The resulting section drag 
coefficients were measured using a wake survey probe. Two 
flight regimes were examined during the test; 1) cruise, 
with high velocity and low angle of attack, and 2) climb, 
with low velocity and high angle of attack. The 
temperature in the tunnel was set to -4 degrees C to 
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generate glaze ice shapes and -26 degrees C for rime 
shapes. 

Two methods are available for recording the ice 
accretion geometry. For short icing times, a small 
section of ice is scraped away near the leading edge of 
the model. A template is then inserted into the gap and a 
tracing can be made. For longer periods of accretion, a 
section of the ice is removed by spraying steam inside the 
model near the leading edge. It is then dipped into a 
container of molten beeswax. After hardening, the water is 
removed, the plaster is poured inside and casts are then 
available for more detailed tracings [13]. 

From the shapes generated during this tunnel entry, 2 
rime and 2 glaze shapes were chosen to represent typical 

clinb and cruise conditions. These shapes were then 
modelled for the second tunnel entry. Table 1 gives a 
summary of the pertinent test parameters which generated 
the chosen shapes. 


TABLE 1 

Ice Generation Test Parameters 


TYPE 

T 

a 


d 

LHC 

t 

p ice 

BIME 

-26 

2.6 

51 

15 

1.5 

15 

0.421 

BIHE 

-26 

6.6 

40 

15 

1.5 

15 

0.534 

GLAZE 

-4 

2.6 

51 

15 

1.5 

15 


GLAZE 

-4 

6.6 


20 

2.9 

15 



A fifth shape, denoted Generic Glaze was derived from 
the work of Ingelman-Sundberg £12], This shape was chosen 
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because it readily scales down tc a 6" chord model. 
Comparison testing of this shape will be performed in the 
Ohio State Transonic Wind Tunnel facility. 

The simulated ice shapes were formed from mahogany and 
extended full span. In order to obtain surface pressures, 
the inside of each shape was hollowed out to allow 
clearance for the 1/8" ID tubing required for tapping 
(Figures 5-7). 

In order to obtain pressures on the airfoil itself, 
1/8" OD strip- a-tube was attached to the surface. In 
order to simulate the natural roughness of ice accretions, 
aluminum oxide grit with a k/c=. 00058 was attached using 
an acrylic spray adhesive to the glaze shapes, while a 
grit with a k/c=.0012 was added to the rime shapes. 

Data acquisition and reduction was accomplished using 
the OSO Digital Data Acquisition and Seduction System 
(DDABS) [20], The heart of the system (illustrated in 
Figure 8) is a DEC ISI-11 microcomputer. System input and 
output is through a standard teletype terminal, and the 
mass storage device is a single-head dual-drive floppy 
disc. Signals from the various pressure transducers and 
the wake probe slidewire enter the analog front-end, which 
conditions the signal and converts it into digital format 
for direct input to the microcomputer. 

A Scanivalve transducer system was used to provide 



surface pressures on the model and a twin-head wake survey 
probe, with wake total and static ports, was used to 
sample pressures in the wake. Drag data were then obtained 
using the wake momentum deficit technique. Figure 9 shows 
a schematic of the data acquisition system set-up. 

One of the key features of the 0S0 DDARS is on-line 
data reduction. The system operator is given quick-look Cp 
distributions as well as integrated values of C^, Cm, and 
Cd. The engineer can then evaluate the progress of the 
test and maximize tunnel usage time. 

Final data reduction was performed cn the OSU Harris/6 
Computer System. Hard-copy plots of the Cp distributions 
for each configuration were generated and integrated 
values of lift, moment, and drag coefficient were 
obtained. 

In order to visualize the flow in the region of the 
ice shape, a splitter plate arrangement was used. The 
plate could be inserted into place between the upper and 
lower halves of the simulated ice shape (Figure 10) . 
Droplets of oil-based paint were applied and the tunnel 
brought from idle up to the required test speed, after no 
further movement of the droplets was observed, the tunnel 
was brought to idle and photographs were taken of the flow 
patterns (Figures 11-13). The separated streamline 
coordinates were digitized from these photographs for use 
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in the mixed analysis/design study. 

Five configurations were run during the two year 
program, including deflecting the flap from 0-30 degrees. 
Of importance to this report were results obtained on: 

1) . Glaze 3 

2) . Glaze 7 

3) . Generic Glaze 

The Cp distributions and integrated lift coefficients 
provided the necessary database for the analysis effort 
which will be discussed in Chapter 5. 
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V. BESUITS AND DISCUSSION 


luc approaches to evaluating the performance of a 
glaze ice shape were used in this study. Ecth relied upon 
the database generated in the lewis Icing Tunnel on the 
simulated ice shapes. The first scheme was tc examine 
current airfoil analysis codes and compare the predicted 
inviscid pressure distribution to the experimental result. 
The second approach utilized the Bristow inviscid design 
and analysis program in an attempt to predict the shape of 
the separated zone behind the glaze ice horn. Together 
with this effort equivalent tody concepts were 
investigated. 

Analysis of Current Potential flow Schemes 

As a first attempt at analyzing glaze ice accretions, 
an investigation of current airfoil analysis programs was 
performed. Computer programs utilized in this phase were 
Smetana, Eppler, Dvorak, Eristcw, and Theodcrsen. 

To initially evaluate these programs, sample cases 
were run on the clean 63A415 airfoil and compared to 
experimental results obtained in the Lewis Icing Tunnel. A 
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representative comparison is shove in Figure 14. This 
particular distribution is at a=2.6° and vas obtained from 
the Bristow code, tut the results cf all the programs 
studied were nearly identical. Good agreement with 
experiment vas seen. 

An interesting observation can be nade about the 
various panelling methods described in Chapter 3 of this 
text. Ihroughout this phase of the study, very little 
difference vas seen among the pressure distributions 
generated by the Eppler and Bristow programs. However, Cp 
distributions from the Eppler analysis do show a higher 
degree of sensitivity to the coordinates. This can be seen 
in the higher freguency and magnitude of pressure spikes, 
particularly in the leading edge region. This is primarily 
due to the means by which Eppler cubic splines the input 
coordinates to define the panels. Smoothing cf all ice 
shapes was a necessity for input tc this program. 

The first ice shape tc be examined, the Glaze 7 case, 
was a logical progression from the clean airfoil. As seen 
from Figure 6, this shape is ncnctcnically increasing in 
X. Figure 15 shows a resulting pressure distribution from 
the covparisons made. Prediction again is very good at 
this angle cf attack, 4.6°. However, as the angle of 
attack was increased and the laminar separation babble in 
the region cf the ice shape herns grew, the potential flow 


results were not very good, This is understandable in 
light cf the highly viscous nature of the separation 
bubble. 

Three of the studied computer programs had boundary 
layer routines; Dvorak, Eppler, and Smetana. However none 
had the capability to predict the separation bubble 
geometry and flow properties. Shen laminar separation was 
predicted, the bubble was assumed to be small enough to be 
considered negligible. Thus re-attachment was predicted at 
the same location as separation. Ibe flow was then 
considered turbulent from this point on. However, due to 
the large adverse gradient in this area the turbulent 
boundary layer routines soon predict separation also. It 
should be noted that the laminar separation point 
predicted by Dvorak compared very well with the observed 
flow visualization separation point. 

Flow visualization techniques however reveal the true 
size of the separation bubble (Figures 11-13). Bubble 
lengths of 10% chord were observed at moderate angles cf 
attack. This definitely shows that the assumptions made by 
these computer programs, even though valid for most cases, 
break down when applied to the flow in the region of the 
ice growth. 

The Glaze 3 and Generic Glaze shapes, due to the fact 
that they are not mcnctcnically increasing in X, proved to 
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be much note difficult tc analyze* The Smetana program 
simply would not run on a double-valued shape and the 
Theodorsen conformal napping method (Figure 16) could not 
successfully map the iced airfoil to the circle plane. 

another difficulty arose at this time with the Dvorak 
program. Figures 17-19 show the panel geometry produced by 
the Bristow, Dvorak, and Eppler cedes respectively, for 
the Glaze 3 case. Hbile Bristcv and Eppler modelled the 
large change in slope very well (Eristow dees not 
redistribute the coordinates) , Cvcrak«s method poorly 
approximated the geometry. Figure 18 shows a panelling 
attempt by Dvorak for the Glaze 3 shape. The lower horn 
was net retained in the panelled configuration. This 
inability to correctly represent the input geometry was 
seen throughout the analysis of the Glaze 3 and Generic 
Glaze shapes. 

Figures 20-22 show the comparisons between theory and 
experiment for the Glaze 3 ice shape at a low angle of 
attack. Reasonable accuracy is obtained for this case. 
However, when the angle of attack was increased, results 
degraded quickly. Figures 23-25 show the Generic Glaze 
shape at a moderate angle of attack, 5.6°. As the angle of 
attack of the airfoil with ice is increased, the viscous 
effects become guickly much more important than for clean 
airfoils at a similar angle. 
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Figure 23 shows the difficulty associated with trying 
to treat a viscous flow problem with an inviscid approach. 
The large pressure spike, observed with all ice shapes, 
occurs at the tip cf the horn as the flow attempts to 
negotiate the large change in surface slope at this pcint. 
None of the programs examined could predict the observed 
constant pressure zone associated with the laminar 
separation bubble. 

Even though comparisons between theory and experiment 
made at low angles of attack were good, when moderate 
angles are evaluated the viscous effects associated with 
the ice shape need to be considered. Table 2 shews this 
very clearly. It should be noted that the theory row 
corresponds to an averaging of the results from Eristcw, 
Dvorak, and Eppler for that angle of attack (Smetana was 
included for the Glaze 7 cases), figure 26 shows a summary 
of the characteristics of the airfoil analysis methods 
investigated. The next step in the analysis then was to 
examine the shape and length of the laminar bubble. 

TABLE 2 

Lift Coefficient Prediction with Ice 
GLAZE 3 

a -2.4 3.6 5.6 9.6 
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1.57 

1.18 


€.6 

1.45 

1.30 


3.6 5.6 

0.85 1.10 

0.72 0.84 

Mixed Analysis/Design Method 

The Bristow program has the unique option of 
performing mixed analysis and design problems. This 
feature was utilized in an effort to predict the shape of 
the laminar separation zone. 

The input to the Bristcv mixed analysis/design option 
involved holding the geometry fixed at the tips of the ice 
horn (1 panel was fixed on each hern). In addition, 
tangential and normal velocities in the design region were 
reguired. all normal velocity components sere set to zero. 
The tangential component was then calculated from the 


Theory 

0.10 

0.84 

1.09 

Experiment 

0.08 

0.75 

1.01 



GIJZI 

7 

a 

-3.4 

2.6 

4.6 

Theory 

-0.03 

0.72 

0.96 

Experiment 

-0.03 

0.70 

0.90 


GESEBIC 

GLAZE 

a 

-2.4 

-0.4 

1.6 

Theory 

0.10 

0.35 

0.60 

Experiment 

0.10 

0.32 

0.54 


29 



experimental pressure coefficients and Eernoulli*s 
Equation. 

Since quantitative flow visualization data was cnly 
availatle for comparison for the opper surface, the 
geometry of the buttle in this region was studied 
primarily. However, conclusions drawn here should apply in 
the lower surface separated zone and the reqion between 
the two ice horns, from the photographs of the splitter 
plate arrangement (Figures 11-13), digitized coordinates 
for these regions were obtained for comparison to theory. 

One final parameter needed to te examined before 
prediction of the bubble geometry could be made. This 
parameter, the reattachment point, is the position on the 
airfoil up to which velocities are specified and beyond 
which geometry is fixed. Figure 27 shows the predicted 
geometry of a separation bubble on the Glaze 3 shape. The 
reattachment point was varied from X=.04 to X=.80 . The 
shape of the bubble converged to the solid line in this 
figure. Moving this point further tack on the airfoil 
surface did not alter the shape of the bubble. Therefore, 
for the cases examined here rear position of the design 
region was set to X=.20 . 

figure 28 shows a comparison between predicted and 
experimental shapes of the separation bubble. Seasonable 
agreement is seen at this lew angle of attack. However, as 
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the angle of attack was increased, the predicted shape 
tended to he longer and thicker than the observed one. 
Figure 29 shows a comparison run on the Glaze 3 shape at 
5.6° . Experimentally, the reattachment point was 
observed to be at X=.Q5 . Thecret ically however, it was 
found to be at X=.175 . 

There are a number of reasons for these discrepancies. 
First, with a splitter plate technigue cf this kind, the 
line that is visualized is actually a little above the 
zero velocity line (Figure 3) , not the separated 
streamline. This would agree with the observation that the 
splitter plate shape lies within the bounds of the 
theoretical prediction. 

A second, and far more important difficulty was 
discovered while studying the flow visualization 
photographs. In Figure 11, the streamlines are observed to 
converge, indicative of a flow no longer 2-D in nature. A 
test program was performed in the OSU Subsonic Sind Tunnel 
to determine the nature of this problem [21]. 

A GAS-1 airfoil was outfitted with a splitter plate 
and a simulated glaze ice shape. The airfoil was run 
through a series of angles cf attack, first with the 
splitter plate leading edge protruding out into the 
stream, and second with this portion of the plate removed. 
The results of this study show that with the larger 
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splitter plate, the boundary layer separates off the plate 
and induces vortices due to the impressed adverse pressure 
gradient from the ice shape. These vortices traveled 
downstream, affecting the 2-D nature of the flow near the 
splitter plate. Quantitative measurements showed a change 
in re-attachment point cf 5* was possible between the two 
plates. This value however cannot be directly applied to 
the results on the 634415 airfoil in the IBT. Bather, the 
reader should use this information qualitatively when 
applying it to Figures 28-29. The important point is that 
the large splitter plate moved the reattachment point 
forward on the airfoil surface. Keeping this factor in 
mind the prediction of the separated zone in Figures 28 
and 29 appear to be better that first thought. 

A third difficulty with this type of mived-mode 
analysis and design cooes from the assumption that the 
pressure gradient through the boundary layer is 
negligible. This assumption is a bey element of the design 
process but may net be a valid one for the thick 
separation zones associated with glaze ice. 

lastly, comparison with flow visualization is not 
possible in the region between the glaze herns due to the 
reasons just mentioned. However, Figure 30 shows the 
predicted geometry using this method for the Glaze 3 shape 
at a =5.6 °. 
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Egoivalent Body Approach 


The last phase of this study looked at the equivalent 

body approach in which pressures were calculated on the 

input observed separated streamline, Figures 31-33 show 

the pressure distribution in the separated zone for the 

o 

Glaze 3 airfoil at a=5.6 . Bespectively these results are 

from the Bristow, Dvorak, and Ippler codes. The dashed 
lines represent an inviscid solution obtained from the 
physical airfoil geometry only. The solid lines are the 
iaproveaent obtained when the coordinates of the separated 
streamline from the flow visualization are input. The 
improvement does net appear very significant for this case 
but that is primarily due to the position and extent of 
the bubble. It should be noted that the coordinates of the 
separated streamline were not smoothed before input. As a 
result, a large pressure gradient is obtained where the 
separated streamline rejoins the airfoil surface (Figure 
29) . 

Figures 34-35 show another comparison with a thicker 
and longer separation bubble. With this case, a vast 
improvement is obtained between the inviscid prediction 
based on the actual geometry and that based on the 
separated streamline. Particular notice should be taken of 
the comparison in the area of the separated zone behind 
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the ice shape hern. Lastly, a test was performed of the 
design method of the Bristow code using these conditions. 
The pressure distribution calculated by Bristow for the 
separated streamline. Figure 34, was re-input as a design 
region. The geometry predicted from this distribution is 
shown in Figure 3€ along with the original separated 
streamline geometry. Excellent agreement is obtained and 
substantiates the use of the Bristcw program for these 
applications. 
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VI. SUHHABI ABE COBCltJSIOHS 


An experimental program was conducted to expand the 
current datatase of performance data on airfoils with 
glaze ice. simulated ice shapes were developed based on 
actual ice growth on the BACA 63Ai»15 airfoil in the HASA 
Lewis Icing Besearch Tunnel. These shapes were tapped so 
pressure distributions could be obtained. In addition, 
flow visualization photographs were taken of a splitter 
plats arrangement in the region around the ice shape. 

Extensive comparisons were run using current airfoil 
analysis programs such as Eppler and Smetana in an effort 
to predict the pressure distribution and separation zone 
geometry of these ice shapes. Also, comparisons were made 
using the Bristcw Hixed Analysis/Eesign Program between 
the separated streamline geometry obtained from the flow 
visualization and the predicted geometry designed from 
input values of velocity. The following conclusions can be 
made from the study described here: 


1. Host panelling methods can predict the pressure 
distribution of an airfoil with ice, tut only at a low 
angle of attack. Bhen the angle is increased to 
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noderate levels, the method breaks down because of the 
large separation bubble created and its viscous 
nature. 

2. Panelling methods that do have boundary layer 
routines treat the laminar bubble as a transition 
point from laminar to turbulent flov. This transition 
is considered to occur in a negligible distance. 

3. The classical assumption that the pressure gradient 
through the boundary layer is negligible appears tc 
bold even for the thick separation senes associated 
with glaze ice accretions. Seasonable predictions of 
the bubble length and shape were obtained from this 
assumption. 

4. Improved results are obtained from the theory when 
an equivalent body approach is applied. The 
coordinates of the separated streamline are input 
rather than the physical geometry of the airfoil 
surface. 

It is recommended that before an attempt is made tc 
develop a numerical approach tc analyze glaze ice 
accretions, the following steps are taken: 

1. Obtain more detailed pressure distributions in the 
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separated zone behind the ice horns and between then. 
The more detailed the surface pressure distribution 
is, the better the results the nixed analysis and 
design progtan yields. 

2. Obtain pressures vertically through the separated 
zone, also, at the sane tine neasure the velocity 
profile in this region. lastly, a deternination should 
he aade cf the transition point from laminar to 
turbulent flow in the shear layer. 

3. Repeat the splitter plate flow visualization 
experiments with a snaller plate so as not to ruin the 
2-D nature of the flow. Ihis will give a tetter idea 
where the reattachnent point is. 
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OSEF’S GUIDE TO THE E5IST01 CODE 


This chapter is intended as a user's guiae for the 
Bristow program. Two modes of operation are possible with 
this program; 1) Analysis only and 2) Bixed Analysis and 
Design. Where applicable the differences in input 
parameters between these modes will be pointed out. 


CARD 1 COLUHNS 1-72 ATITLE 

Enter case title on this card. 

CABD 2 COLUHNS 1-10 ISAfE 

Set ISA7E=0 to indicate the start of a new set of 

2 € <i- isave= 2 for input of a new ALPHA only. Submit 

Set ISAVE= 1 if only retaining QT and N (Q) from 
previous case. All ether inputs can be changed. 
Untransformea (XB,IB) coordinates are reused. 

Set IS A1E=3 to repeat last case wxth new values of 
ALPS A, CIBCE § and VBP distribution. Ho design cases 
are allowed. Only submit cards 1# 2, 6, 8, and 11. 
COLUMNS 11-20 ALPHA 

Angle of attack of x-axis with respect to free stream 
velocity 


CABD 3 COLUHNS 1-10 CT 

Number of airfoil elements (Normally set CT-1. 
flap present set CT=2, etc.) 


If 


CABD 4 COLUMNS 1-10 CBGBD 

Beference length for moment and 
integration (Normally set=1) 
COLUMNS 11-20 CAP P A 1 
Recommend set CAPPA1— .01 . Used 
sharp corner ccntrcl point 
COLUMNS 21-30 CAEEA2 
Becommend set CAEEA2— .02 • Used 


lift coefficient 
in calculation of 
in calculation of 
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Kutta condition control points 
COLUMNS 31-4C LINSIG 
Singularity choices: 

Set LINSIG=0 for constant source distribution on 
panels. 

Set LINSIG= 1 for VINE portion of source distritution 
to be piecewise constant and yNP portion to be 
piecewise linear. 

Set LINSIG=2 for linear source distritution on panels 
(NOTE: LINSIG Choice has little effect on results. 
Becommend set IIHSIG=0 or 2. 

COLUMNS 41-50 VISE 

Non-dimensional free stream velocity (Ncraally set=1) 
COLUMNS 51-6C VBEE 

Ncn-di«ensional reference velocity used to calculate 
pressure coefficients (Normally ?BEf=VIHE) 

CABD 5 COLUMNS 1-10 UMAX 

Number of iterations in design node (Set=0 in 
analysis mode). Suggest set=4 for design node. Most 
cases converge in this number of iterations. 

COLUMNS 11-20 ITB 
COLUMNS 21-30 BIX 

Beconmend set BLX=1.0 . Cesign region geometry is 
relaxed by a factor of ELX every ITB iterations. 
COLUMNS 31-40 ITHICK 

Normally set=0. Allows no thickness increase if 
design process results in negative thickness. 
Execution will terminate if this occurs. Set=1 and 
thickness increase will te allowed inspite of 
negative thickness occurring. 


The following cards should te input for each of the 
elements <Q=1 # C=2, ...» C=CT) 


CABD 6 COLUMNS 1-10 N 

Number of coordinates defining element C 
COLUMNS 11-20 PT 

Number of boundary conditions in element C (PT = 
Number of analysis regions ♦ number of design 
regions) 

COLUMNS 21-30 KUTTA 

Normally set=5. Set=1 to input desired circulation 
normalized ty perimeter of this element. Set=2 to 
input circulation (not normalized ty perimetez). 
Set=3 if Kutta condition is zero velocity normal to 
trailing edge tisector. Set=4 to determine 


39 



circulation from input tangential velocities- Set=5 
for same condition as 3 tut higher order 
extrapolation for trailing edge bisector is used (4 
panels - 2 upper surface, 2 lower surface) 

COLUMNS 3 1 — *1 0 CIBCE 

Input circulation of this elewent. Set=0 for KU1TA>2 
COLUMNS 41-50 DXTE 
COLUMNS 51-60 DUE 

Trailing edge opening. Ignored if trailing edge 
regions are ETIPE=0 (DXTE = X N - X^) (Dili = ! N - 

h) 

CARD 7 COLOURS 1-10 AC 
COLUMNS 11-20 BC 
COLUMNS 21-30 ALEQ 
CCLOMNS 31-40 SCIQ 

Normally set AQ=0, BC=0, ALFC=1# and SC1Q=1. These 
axe transf crmation parameters which are applied to 
input coordinates (XE,1B) to produce a new series of 
coordinates for use in the program. This allows 
translation, rotation, and stretching of the input 
coordinates. The transforaation applied is: 
X-AQ+SCLC*[ XE*COS (ALEQ) -TE*SIS (ALFQ) ] 

I=BQ*SCLQ*[ IB*COS (AIFQ) ♦XE*SIH (ALPQ) ] 

CCLUMNS 41-50 ICLR 

Set=0 for internal flow (counter-clockwise coordinate 
input). Set=1 for external flow (clockwise coordinate 
input) • 


The following cards are input for each of the regions P=1, 
P=2, ..., P=ET of element Q. 


CARD 8 COLUMNS 1-10 MEAN 

Number of panels in this boundary condition region 
(NOTE: KEA*=li-1) 

CCLUMNS 11-20 ISHP 

Set=0 if first point in the region is not a sharp 
corner point 

Set=1 if first point in the region is a sharp corner 
point (NOTE: a sharp corner point is defined as a 
point cf slope discontinuity) 

COLUMNS 21-30 PTIEE 

Set=0 if analysis region with no translation 
Set=1 if analysis region with translation allowed 
Set=2 if design region with first coordinate fixed 
and previous boundary condition was a design region 
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Set=3 if design region with first coordinate free or 
previous boundary condition region was analysis with 
no translation 
CCLDHHS 31-40 PDSf 

Normally set=0. Set=1 for this region to undergo same 
relative length change as previous region. (NCIE:Tbis 
region must he PTXEE=3 and previous regior must he 
P1XPE>2) 

CCLDHHS 41-50 I7HP . . ... 

Set=0 if normal velocities are presented in tbiw 

reai.cn 

Set=1 if normal velocities are all to te set tc zero. 
Ibis is normally the case. 

CCLDHHS 51-60 8 ...... 

Normally set=.001. If large number is input length 

variation is suppressed. 

CABD 9 COLDHHS 1-10 XE 

CCLDHHS 11-20 TP . _ . 

Coordinates of first point in this boundary condition 

region. Ignored if ETXPE£2. 

CCLDHHS 21-30 XBE 
CCLDHHS 31-40 YBB 

If this is a design region and is followed by an 
analysis region, these coordinates are considered to 
be the last point in this region. 

CABD 10 CCLDHHS 1-10 7NE of Panel # 1 

CCLDMNS 51-60° VHP of Panel NPAB 

Omit this card if IVHP=1. Otherwise enter panel 

midpoint normal velocities. 

CCLDHHS 61-70 HBD „ . . . _ , 

Humber of values of VNE on this card. Omit if 6 
values of VHE are on this card or it is the last card 
fer this region. 


HEX1 CABD SEBIES 

CCLDHHS 1-10 

9 

CCLDHNS 51-60 
Omit this card 
panel midpoint 
CCLDHNS 61-70 
Same as HBD cf 


VIP of Panel 1 
♦ 

VIP of Panel HEAR 

if analysis region. Otherwise enter 
tangential velocities. 

NBD 

card 10. 


NEX1 CABD SEBIES 

CCLDHHS 1-10 V TEE of Panel 1 

w • 

CCLDHNS 51-6G VIEE Of Panel NE*N 
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Omit this card if analysis region. Otherwise enter 
panel endpoint tangential velocities* 

COLUMNS 61-70 NBD 
Same as NBC cf card 10 


Input the next series of cards for eleient Q-1, C-2, 
...» C=Q1. emit these cards if IS8VE=0- 


NEXT CfiBD SEBIES . 

cciohns 1-10 XB cf point 1 on this element 

t * 

„ P 


columns 51-60 XE cf point H on this element 
x-coordinates cf airfoil geometry. If external flew, 
input should he clockwise. If internal flow input is 
counter-clockwise. 

COLUMNS 61-70 NED 
Sane as NBC of card 10 


NEXT CftBD SEBIES _ . . 

CCLOWNS 1-10 YB cf point 1 on this €le®€Ct 


COLUMNS 51-60 YB of point N on this element 
I-coordinates of airfoil geometry. If external flow, 
input should fce clockwise. If internal flow input is 
counter-clockwise. 

COLUMNS 61-70 HBD 
Same as NED of card 10 
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PHYSICAL AIRFOIL SHAPE 


VORTEX FILAMENT 



REPRESENTATION OF AIRFOIL BY POLYGON 
AND POINT VORTICES 


FIGURE 2. SMETANA PANELLING METHOD 
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3. DIAGRAM OF LAMINAR SEPARATION BOBBIE 




FIGURE 4. THE NASA LEWIS 6' x 9’ ICING RESEARCH TUNNEL 
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5. GLAZE 3 SIMULATED ICE ACCRETION 
AND PRESSURE TAP LOCATIONS 
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FIGURE 8. OSU DIGITAL DATA ACQUISITION 
AND REDUCTION SYSTEM 
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FIGURE 9. OSU DATA ACQUISITION SYSTEM 
AS USED IN THE NASA LEWIS IRT 
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FIGURE 10. SPLITTER PLATE ON 63A415 AIRFOIL 
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Upper Surface 
M= 0.152 , 

Re = 4. 7 X 10 b 



FIGURE 11. SPLITTER PLATE PHOTOGRAPH CF GENERIC GLAZE 
ICE SHAPE ON 63A415 AIRFOIL (o=-0.4°) 




Upper Surface 
M= 0.152 
Re = 4. 7 X 10 6 




FIGURE 13. SPLITTER PLATE PHOTOGRAPH OF GLAZE 3 ICE 
SHAPE ON 63A415 AIRFOIL (a=5.6°) 
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IRT Test (C p =. 85) 
a=4. 6° * 

M= . 123 c. 

Re=4. 0 x 10 

Bristow (C£=.96) 
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FIGURE 15. COMPARISON BETWEEN EXPERIMENT AND 
THEORY FOR THE 63A415 AIRFOIL WITH GLAZE 7 

ICE SHAPE 
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FIGURE l€. TIIEODORSEN TRANSFORMATION OF 63A415 
AIRFOIL WITH GLAZE 3 ICE SHAPE 
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FIGURE 18. DVORAK PANELLING SCHEME FOR 
GLAZE 3 ICE SHAPE 
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FIGURE 19. EPPLER PANELLING SCHEME FOR 
GLAZE 3 ICE SHAPE 
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FIGURE 24. COMPARISON BETWEEN EXPERIMENT AND THEORY 
FOR THE 63A415 AIRFOIL WITH GENERIC GLAZE ICE SHAPE 
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CODE 


POTENTIAL SOLUTION 


RUN TIME 


COMMENTS 



EPPLER 

MIXED PANEL METHOD 
PARABOLIC VORTICITY 

2 MIN. 

EXTREMELY SENSITIVE TO 
GEOMETRY 

REQUIRES ICE SHAPE 
SMOOTHING 

SPLINE FITS TO FORM 
PANELS, C f METHOD 

MAX 


SMETANA 

PANEL METHOD 
CONSTANT VORTICITY 

2 MIN. 

X MONOTONICALLY 
INCREASING 


DVORAK 

PANEL METHOD 
LINEAR VORTICITY 

2 MIN. 

REDISTRIBUTES AIRFOIL 
COORDINATES 

POOR ICE SHAPE MODELLING 
C p METHOD 

MAX 


BRISTOW 

PANEL METHOD 
SOURCE AND VORTICITY 

5 MIN. 

RELATIVELY INSENSITIVE 
TO ICE GEOMETRY 
MULTI-ELEMENT MODE 
DESIGN WITH MIXED BC 


WOAN 

THEODORSEN 
CONFORMAL MAPPING 

1 MIN. 

SENSITIVE TO GEOMETRY 


FIGURE 26. FLOWFIELD PREDICTION METHODS SUMMARY 





Reattachment point®. 08 
Reattachment point®. 15 
Reattachment point®. 20 



FIGURE 27. SEPARATION ZONE PREDICTION FROM MEASURED Cp's FOR THE 
63A415 AIRFOIL WITH GLAZE 3 ICE SHAPE AND VARYING REATTACHMENT POINT 



Calculated Separation Zone from 
Measured Cp's Using Bristow Code 


Measured Separation Zone from 
IRT Flow Visualization 



FIGURE 28. SEPARATION ZONE PREDICTION FROM MEASURED Cp’s 
FOR THE 63A415 AIRFOIL WITH GENERIC GLAZE ICE SHAPE 



Calculated Separation Zone from 
'Measured Cp’s Using Bristow Code 



FIGURE 29. SEPARATION ZONE PREDICTION FROM MEASURED Cp's 
FOR THE 63A415 AIRFOIL WITH GLAZE 3 ICE SHAPE 
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FIGURE 30. PREDICTION OF REGION BETWEEN GLAZE ICE HORNS FROM MEASURED Cp ' s 
FOR THE 63A415 AIRFOIL WITH GLAZE 3 ICE SHAPE 





FIGURE 32. PRESSURE DISTRIBUTION IN SEPARATED ZONE 


BEHIND UPPER SURFACE HORN OF GLAZE 3 ICE SHAPE 
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FIGURE 33. PRESSURE DISTRIBUTION IN SEPARATED ZONE 
BEHIND UPPER SURFACE HORN OF GLAZE 3 ICE SHAPE 
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IRT Test a=3 . 6 
M= . 154 , 

Re=4.8 x 10 b 


Potential Flow Result 
Based on Actual Airfoil 
Geometry (Bristow) 



FIGURE 34. PRESSURE DISTRIBUTION IN SEPARATED ZONE 
BEHIND UPPER SURFACE HORN OF GENERIC ICE SHAPE 


IRT Test ot=3.6 v ' 

M=.154 , 

Re=4 . 8 x 10 b 

Potential Flow Result 
Based on Actual Airfoil 
Geometry (Dvorak) 

Potential Flow Result 
Based on Separated 
Streamline from Flow 
Visualization (Dvorak) 


FIGURE 35. PRESSURE DISTRIBUTION IN SEPARATED ZONE 
BEHIND UPPER SURFACE HORN OF GENERIC ICE SHAPE 




Predicted Separation Zone from 

Cp's Calculated for the Flow Visualization 

Shape 

— Measured Separation Zone from 

IRT Flow Visualization 



FIGURE 36. EVALUATION OF BRISTOW DESIGN METHOD 
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